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Realistic Surface Scattering and Surface Bound State Formation in the High T c 

Superconductor YE^CusOg+x 
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Surface Umklapp scattering of quasiparticles, and surface roughness are shown to play essential 
roles in the formation of the surface bound states in realistic models for YBa,2Cu^O(,+ x . The results 
account for the shape, the impurity dependence of the height, and for a proposed universal width 
of the zero bias conductance peak. 

PACS numbers: 74.20.-z, 74.25.Jb, 74.50. +r, 74.80.fp 
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This article describes the formation of mid-gap sur- 
face bound states in models of YBa 2 Cu306+2; (YBCO) 
which include realistic surface scattering of the quasipar- 
ticles. (We define realistic surface scattering to be scat- 
tering that includes surface roughness as well as surface 
Umklapp scattering processes and a realistic Fermi sur- 
face geometry.) More than a trivial extension of what is 
known about surface bound states from the study of el- 
ementary models is required. (By our definition, 
elementary models neglect either surface roughness or 
surface Umklapp scattering or both). On the one hand, 
except at grazing angles, the superconducting quasipar- 
ticles are prevented from forming surface bound states 
because they are reflected diffusely at a rough surface, 
and on the other hand, the grazing incidence quasipar- 
ticles are prevented from forming surface bound states 
on a [110] surface (the most favorable case) when sur- 
face Umklapp scattering is considered || . Surface bound 
states are of interest in the study of high T c superconduc- 
tors, because they probe the sign changes of the gap func- 
tion A(k) on the Fermi surface, and also because they 
manifest themselves experimentally through the zero bias 
conductance peak (ZBCP) observed in tunneling exper- 
iments. As noted in the conclusions, our results have a 
considerable impact on the interpretation of experiments. 

The existence JlJ of surface bound states in d-wave 
superconductors and their role in producing the ZBCP 
observed in tunneling experiments fi"(i|-|l8|| are well estab- 
lished for elementary models, e.g. in g]. Further stud- 
ies have discussed, among other results, the existence of 
multiple sub-gap resonances || , the surface-induction of 
order parameters of different symmetries [j^]7|,|il| and 
the associated splitting of the ZBCP at low tempera- 
Q, the splitting of the ZBCP in a magnetic 



13 , the effects of irradiation in suppressing the 
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According to the Rayleigh criterion in the theory of 
surface roughness (for a general discussion see |2l[] ) , a sur- 
face is rough (reflects diffusely) or smooth (reflects spec- 



16], the effects of surface roughness in broaden- 
and the anisotropy of the ZBCP 



ularly) according as the parameter Rk = k±ij is smaller 
than or greater than unity. Here k± is the component 
of the incident-wave wave vector perpendicular to the 
surface and rj is the surface asperity. For Rk small, per- 
turbation theory (e.g. see pi| ) gives the probability of 
diffuse reflection as Pd ~ i.e. it goes to zero as k\ 
for small k±. On the other hand, for Rk large, a quasi- 
classical approximation shows that, for a wave incident 
normally on the surface, the probability of specular re- 
flection is given approximately by P$ = exp(—4irRl), i.e. 
there is essentially no specular reflection for Rk greater 
than unity (e.g. see p^]). 

It follows from the above estimates that, except for 
quasiparticles approaching the surface at grazing inci- 
dence, the reflection is totally diffuse. These ideas have 
been successfully applied to the theory of the electronic 
surface states of a normal metal in low magnetic fields 
p3| where microwave absorption experiments show that 
quantized magnetic surface levels exist only for electrons 
incident on the surface at grazing angles (so that they are 
specularly reflected). The application of these ideas to 
surface bound states in d-wave superconductors suggests 
that, except possibly for those which approach the sur- 
face at grazing angles, the quasiparticles will be diffusely 
reflected at a rough surface and thus will not maintain 
the phase coherence necessary (e.g. see the argument 
surrounding Eq. 25 of to form bound states. 

Surface bound states in a model which incorporates 
surface Umklapp scattering have recently been described 
in 0k As in |||4j, to obtain a realistic Fermi surface 
|p5|^6| we use a discrete lattice model (see Fig. |l|) in- 
cluding both nearest and next nearest neighbor hopping 
interactions in the superconductor. Also, we consider the 
case of a [110] surface, for which the formation of surface 
bound states is the most favorable. The surface bound 
states are eigenstates of the two-component Bogoliubov- 
de Gennes equations and are linear combinations of basic 
solutions of the form 



Ck,(x) = 



A, 
E±iQ,i 



(i) 



l 



as in HJHI^ and others. Here ki = (k ix ,k y ) must be 
on the Fermi surface, Aj = A(ki), = yAf — E 2 , and 
Ki = Qi/(h\vi x \) with Vi X being the x component of the 
normal-state electron velocity. The upper and lower signs 
in Eq. [j] correspond to vi x < and Vi X > 0, respectively. 
From Fig. || it is seen that for k y > k y , a line of constant 
k y intersects the Fermi surface at four distinct points, 
giving four values ki x , i— 1 ... 4, and hence four lin- 
early independent exponentially decaying solutions. In 
this case the surface bound state solutions have the form 



*7(x) =£ i C i £/ ki (x) 



(2) 



From above, only grazing incidence quasiparticles, i.e. 
those corresponding to small values of \ki x \ in the surface- 
adapted Brillouin zone || of Fig. |^, will have the possi- 
bility to form surface bound states in the case of rough 
surfaces. However, it turns out that a detailed study of 
such states |jj has shown that, at least in the case of per- 
fectly fiat surfaces, waves having k y > ky in Fig. || (and 
this includes all waves of small can not form sur- 

face bound states [the reason has to do with the relative 
signs of the gaps A(ki x , k y ) and A(fc 2 a;, k y )\. 

In summary, on the one hand, except at grazing inci- 
dence, quasiparticles (or holes) striking the surface are re- 
flected diffusely and hence do not have the coherence nec- 
essary to form surface bound states, while on the other 
hand, grazing incidence quasiparticles are prevented from 
forming bound states in a realistic model for the Fermi 
surface because they do not satisfy the necessary condi- 
tions on the sign variation of the gap A(k). 

We now give a qualitative explanation of the result 
of our quantitative calculation presented below showing 
how it is possible to form surface bound states in mod- 
els of YBCO which include realistic surface scattering. 
First, assume that the surface at x = is perfectly flat, 
and consider an incoming wave of the form of Eq. |l] with 
(the x- component of its) wave vector equal to —k± x , 
or equivalently —qi, of Fig. ^. In order to satisfy the 
boundary conditions at the surface, two outgoing waves, 
a specularly scattered wave at k\ x and an Umklapp scat- 
tered wave at k 2x , are required. (Note that the wave at 
k 2x appears as a normal scattering process when viewed 
in the surface- adapted Brillouin zone of Fig. ||, but is a 
surface Umklapp process when viewed in the normal bulk 
Brillouin zone ||.) It is the Umklapp scattering into the 
state with wave vector k 2x which prevents the formation 
of the surface bound state because the relative signs of 
the gaps at wave vectors ±ki x and ik 2x are incompati- 
ble with the formation of a surface bound state |J . Thus 
we can only get surface bound states (or at least reso- 
nances) if the scattering from wave vectors ±k\ x to ±k 2x 
is sufficiently reduced. 

It is the roughness of the surface which limits the Umk- 
lapp scattering of quasiparticles of wave vectors ±k\ x into 
the states of wave vector ±k 2x . Quasiparticles in states 



which are linear combinations of states with wave vec- 
tors ±k 2x will be strongly diffusely reflected by a rough 
surface and will thus have short lifetimes. Their short 
lifetimes will spread out their density of states in en- 
ergy, reducing the density of states at the energy of the 
quasiparticles of wave vector k\ x and hence reducing the 
probability of scattering from wave vectors ±fci x to ±k 2x . 

We now develop a more quantitative theory of the 
effects of surface roughness by assuming that surface- 
bound-state component states (Eq. |l|) having wave vec- 
tors iki X have a finite lifetime Tj = Ti/"fi due to nonspec- 
ular scattering at the surface. For quasiparticles with 
vectors ik 2x (which have a probability of unity for dif- 
fuse scattering at the surface) we take this lifetime to 
be the time taken for an electron to get to the sur- 
face from a distance of (k.;) -1 , which, for \E\ <C IA2I 
gives 72 = |Aa|. (Similar qualitative estimates of damp- 
ing rates were shown in Q to be in excellent agreement 
with the results obtained from a detailed solution of the 
Bogoliubov-de Gennes equations.) For grazing incidence 
quasiparticles of wave vector ±ki x , which have a very 
small probability Pi of diffuse scattering at the surface, 
the appropriate value of 7 is 71 = Pi|Ai|. (The perturba- 
tion theory estimate of Pi is given in the above general 
discussion of surface roughness.) These lifetime effects 
are incorporated into our model by using the phenomeno- 
logical approach ( |7],[l4|] ) of replacing E by E + iji in Eq. 

1 for £/ki . If we now look for solutions of the form of Eq. 

2 which satisfy the appropriate boundary conditions for 
a superconductor to vacuum surface (corresponding to 
the absence of the normal metal and the insulating layer 
in Fig. [l]), and which are valid in the limit of small \ki x \, 
we find a solution having an energy E = Eb + iF r where 
the bound state energy is Eb = and the width due to 
roughness is Tr = 71 + I12, with the component 

T 12 (k y ) = |Ai|(T 2 + \A 2 \)\sin qi cot(q 2 /2)\/y 2 . (3) 



Here T 2 = \J A\ + j 2 and q, L = ki X aj\J~2. Note that Ti2 
depends on k y through the dependences of qi and Aj on 
k y . From the above result Tr = 71 + Ti 2 , we interpret 
the solution as being primarily a state formed from the 
grazing incidence wave vectors iki x and having a width 
71 which is augmented by Fi 2 due to transitions from 
wave vectors ±ki x to ±k 2x . It is clear that for 72 having 
the value |A 2 | suggested above, the width T 12 is smaller 
than j Ax I (since qi is small) and the bound state thus has 
a spread of energies lying within the gap |Ai|. If 72 were 
to be small on the other hand, Ti2 would be large and 
no well defined state would exist in the gap. Clearly it is 
the broadening of the quasiparticles at wave vectors ±k 2x 
by diffuse scattering that allows the formation of surface 
bound states for the case of realistic surface scattering. 

We now proceed to a calculation of the tunneling con- 
ductance for the model NIS junction illustrated in Fig. |l| 
using a procedure well-established in other work, i.e. we 
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use the BTK formalism |27] but extended to the case of 
d-wave superconductors [H|J§] an d to a discrete-lattice 
model P,pO|, The procedure is straightforward: assume 
a solution in the form of a linear combination of an elec- 
tron moving towards the junction together with a re- 
flected electron and reflected hole in normal state, and a 
linear combination of the form of Eq. ^ in the supercon- 
ducting state. The coefficients in the linear combinations 
are found using the boundary condition implicit in Fig. 
|l| and the tunneling conductance is calculated following 
BTK p|,p7||. We restrict ourselves to the low tempera- 
ture limit, and to the calculation of the contribution of 
the surface bound states to tunneling conductance, which 
is valid only for voltages such that eV is less than the 
maximum gap, and we also work in the weak transmis- 
sion limit where the probability of an electron tunneling 
across the barrier is much less than unity, since it is only 
in this limit that one obtains a sharply defined ZBCP in 
the gap |,§. In this way, we find that the conductance 
per surface unit cell (for eV smaller than the maximum 
gap) is given by the formula 



r- 4e 7r (k ) TT{k ^ 



(4) 



where the angular brackets indicate an average over k y . 
Here 

T T (k y ) = l\{k y ) +y tmp + ^\2{k y ) + TsN(ky) (5) 

where the damping rate of the basics states with wave 
vectors ±ki x , called 71 above, has been replaced by 71 + 
limp, thus separating the momentum-dependent diffuse 
surface scattering part 71 (proportional to qf ) from the 
momentum- independent bulk impurity scattering part, 

limp ■ Also 

T S n = 2a N acsinq(sinq 2 jiT 2 - smgi|Ai|7 2 )/72, (6) 

with ctN — tjy /Vb, Qc = tc/Vo, Vq being the additional 
potential on the ions in the insulating layer I of Fig. 
and q > being the wave vector of the incident normal 
state electron. Our calculations are valid to the lowest 
nontrivial order in ajy, ac (i.e. a weak transmission insu- 
lating barrier) which means that the contribution of T$n 
to Tt should be neglected (but that T$n should remain 
in the numerator). The fact that T$n is nevertheless one 
of the contributions to Tt allows its interpretation as a 
contribution to the lifetime of the surface bound state re- 
sulting from the tunneling of the bound state excitation 
from the superconductor into the normal metal. 

The tunneling conductance described by Eq. [| is plot- 
ted in Fig. U as a function of eV/Ag (charge times voltage 
over maximum gap). The contribution of 71 to Tr is ne- 
glected since it is smaller than Ti2 at small wave vectors 
qi . Thus the height, width and shape of the curve of con- 
ductance versus voltage are determined by Ti2, T$n, and 



limp- In particular, the cusp-shaped maximum at eV—0 
in the curve for ii rnp = is due to the fact that Fi2 and 
Tsn both go to zero as q\ at small q%. This means that 
the bound states with very small q\ are very narrow and 
contribute all of their weight near eV — 0, giving a sharp 
maximum there. The broader bound states at larger val- 
ues of qi give broader contributions to the conductance 
and are responsible for the weight in the wings. A cusp 
shaped ZBCP has been observed in | Q . 

A strong reduction of the ZBCP height with increasing 
disorder (i.e. increasing limp) without, however, a signif- 
icant increase in the ZBCP width has been found in [|l6| . 
This agrees with the behavior which is predicted by our 
Eq. |] and shown in Fig. |^. The qualitative explanation 
for this can be found in the momentum dependence of 
Ti2. The very small q\ quasiparticles, which are very 
narrow in energy, give a high ZBCP in the absence of 
impurity scattering. Impurity scattering broadens these 
quasiparticles and reduces their contribution to the peak 
height, but has little effect on the broader quasiparticles 
at less small qi which are responsible for the overall width 
of the ZBCP. 

Eq. (4) predicts a universal width, independent of sur- 
face asperity and impurity scattering (for impurity scat- 
tering not too strong), for the ZBCP assuming that Tt is 
dominated by T^. From Fig. [5[ this universal half width 
at half maximum is approximately Ao /4. This is very dif- 
ferent from the predictions of elementary models where 
the width is due to and increases in proportion to the 
surface asperity or impurity scattering ■ Our 
result is roughly in agreement with experiment where the 
widths of measurements made under widely differing con- 
ditions (for an asperity r\ ~ 1 nm in |Q, for r\ ~ 1 to 10 
nm in Q for rj ~ 500 nm. in |l8[ ] , and for varying ii mp in 
fl6| ) all give aproximately the same half width of about 
2 to 2.5 meV, which is however about a factor of two 
smaller than our prediction. The error in the numerical 
factor is perhaps due to our simplistic representation of 
the insulating barrier by a single atomic layer, or to our 
simplistic treatment of the surface roughness (carried out 
simply by the replacement of E by E + iii). 

Our results, in which the width of the ZBCP is due 
to Ti2, also resolve the problem that fitting the ZBCP to 
elementary models requires an inexplicably large value of 
the so-called smearing factor [see 0(a)]. 

To conclude, we note that the description of surface 
bound state formation and the consequent effects on 
the tunneling conductance in YBCO requires a realis- 
tic description of the surface scattering processes. The 
principal source of width of the surface bound states is 
their limited lifetime due to the surface Umklapp scat- 
tering from grazing incidence states to larger perpendic- 
ular momentum states, and this process hence has a de- 
termining effect on the height, width and shape of the 
ZBCP. These ideas give a qualitative explanation of the 
cusp shaped ZBCP observed in Q, the height reduc- 
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tion (without width increase) of the ZBCP caused by 
increased disorder seen in fj"q| , and the approximately 
universal value (independent of surface roughness and 
impurity concentration) of the width of the ZBCP ob- 
served in IQ^Jpg. 
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done) and the support of the Natural Sciences and Engi- 
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FIG. 1. A model (110) NIS tunnel junction. The solid 
(open) circles, and the solid squares represent ions in the su- 
perconductor (insulator), and the normal metal, respectively. 
Also, tjv is the nn normal metal (and normal to insulating) 
hopping interaction, while ts and tc are the nn and nnn 
superconductor hopping interactions. We also assume a nn 
d-wave gap function, i.e. A a = — A b . 

FIG. 2. The surface-adapted Brillouin zone and Fermi sur- 
face for a (110) surface. The notation q XlV — k x , v a/y/2 is 
used. The Fermi surface shown corresponds qualitatively to 
that measured experimentally (e.g. see ||). Note that q± < 0. 
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